加法比例ハザードモデル

library(survival)

## Loading required package: splines

library(pspline)
# スプラインの自由度の指定した場合
res <- coxph(Surv(time, status) ~ pspline(age, df = 4) + pspline(t5, 
    df = 4), stanford2)
par(mfrow = c(1, 2))
termplot(res, se = T, rug = T)

# AICでモデル選択した場合
res <- coxph(Surv(time, status) ~ pspline(age, df = 0) + pspline(t5, 
    df = 0), stanford2)
par(mfrow = c(1, 2))
termplot(res, se = T, rug = T)

# CAIC(Hurvich et. al. (1996) JRSS-B, 60, 271-293)でモデル選択した場合
res <- coxph(Surv(time, status) ~ pspline(age, df = 0, caic = T) + 
    pspline(t5, df = 0, caic = T), stanford2)
par(mfrow = c(1, 2))
termplot(res, se = T, rug = T)